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ABSTRACT 


A Nusselt-type analysis was performed for laminar film 
condensation on the inside of a rotating truncated cone 
with small half cone angles. This analysis included the 
interfacial shear between the vapor and condensate, the 
vapor pressure drop, the thermal resistance in the condenser 
wall, and the condenser outside cooling mechanism. An 
approximation of the analytical model made it possible to 
find a numerical solution for small half cone angles greater 
than Zero: 

A non-capillary rotating heat pipe containing an 
evaporator, condenser, and distilled water as the working 
fluid was tested. It was rotated at 702 and 1404 RPM, and 
the heat transfer rates of the heat: pipe were determined 
experimentally for different saturation temperatures cor- 
responding to electrical power inputs ranging from 1 kW 
to 9 kW. 

The experimental results showed that the non-capillary 
rotating heat pipe was an effective heat transfer device. 
The approximate numerical solution conservatively predicted 
Cie@iedtatbapsrer rate with a deviation of 18% at 7/02 RPM 


and 5% at 1404 RPM. 
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I. INTRODUCTION 


Pee CONVENTIONAL HERAT PIPE 

A heat pipe is a closed container capable of transfer- 
ring large amounts of heat at nearly isothermal conditions. 
Conventionally it consists of four main parts: (1) evap- 
orator, (2) condenser, (3) working fluid, and (4) capillary 
wicking structure to force the condensate back to the 
evaporator. The working fluid is evaporated in the evap- 
orator section of the heat pipe by the heat source. As the 
vapor pressure increases, the vapor flows toward the con- 
denser end, transporting the heat as the latent heat of 
vaporization. At the condenser end, the heat is removed 
from the vapor as it condenses. The condensate is then 
returned to the evaporator by the capillary pumping action 
within the wicking structure. 

Studies of the heat transfer capabilities have heen made 
on the conventional heat pipe [6, 10]* With the conelusion 
that, for ordinary fluids and moderate heat fluxes, the heat 
pipe is limited in its capability to transfer heat by the 


rate of condensate that can be returned to the evaporator. 


Deol LelatiNG, NON-CAPILEARY HEAT PIPE 
The rate of condensate flow may be increased by removing 


the wicking structure from the condenser. If the condenser 


* . e 
Numbers in brackets indicate references listed in the 
bibliography. i 





end is then internally tapered and the heat pipe is rotated 
about its longitudinal axis, a centrifugal force field is 
set up which will overcome both the viscous losses in the 
system and the gravitational field. 

Ballback's study [10] of the rotating heat pipe indi- 
cated that its heat transfer capability is limited by the 
amount of heat which can be removed from the condenser. In 
his study, Ballback performed a Nusselt-type analysis EOu 
the condensation on the inside of a truncated cone. In 
order to obtain a closed form approximate solution for the 
film thickness, 6, he made the following assumptions: (1) 
o° « tan (large cone angles), (2) no shear stresses between 
the vapor and condensate, and (3) no thermal resistance in 
the condenser wall. Daley [11] was able to improve on the 
theory introduced by Ballback by taking into consideration 
the thermal resistance in the condenser wall and the outside 
cooling mechanism. Furthermore, Daley neglected the assump- 
Eon that gs « tang. Using the velocity profile for the 
condensate developed by allies cise he arrived at a non-linear 
Peco Cader samireren: tal lequation for the film thickness, 
6. Daley numerically integrated this expression using a 
Runge-Kutta-Gill numerical integration scheme with an 
IBM 360 MOD 67 digital computer. To start the integration, 
he used the conditions that 6 = 6i and oo = tangd at X = 0. 
It was postulated that 6i was a function of the over-fall 


condition at the step located at the condenser exit. Daley 


arrived at an expression for the minimum film thickness, 6m, 





at the condenser exit, using the approximation that the flow 
over the step approximated the free fall condition for open 
channel flow. Having found the film thickness from the 
second order differential equation at X = Lc, the condenser 
exit, Daley compared this value with 6m at Le. When the 
difference was less than .004, a solution was found. Using 
this technique, Daley was able to find the heat flow out of 
the condenser as a function of RPM for 0° half cone angle 
only. As the cone angle was increased, however, instability 
of the highly non-linear second order differential equation 
for the condensate film thickness caused the equation to be 
extremely sensitive to the initial value of the film thick- 
A@Ss. Consequently, no solutions were found for the conden- 


sate film thickness for half cone angles greater than zero. 


eels OBIECTIVES 

The objectives of this study were to: (1) continue the 
theoretical study of the heat transfer capability of the 
rotating heat pipe and obtain solutions for $9 > 0°, (2) 
obtain experimental data from the heat pipe designed and 
assembled during Daley's study, and (3) compare the analyti- 


cal and experimental results. 


10 








II. THEORETICAL PROGRAM 


In order to improve on the analysis performed by 
Ballback and Daley the effects of vapor pressure drop and 
interfacial shear stress between the liquid and the vapor 
were included. Furthermore, the thermal resistance in the 
condenser wall and the outside cooling mechanism were con- 
Sidered. Using the coordinate system in Figure 1, for 
small half cone angles, momentum, continuity, and energy 
equations were developed for an infinitesimal fluid element. 
In addition, assuming one dimensional turbulent flow in the 


vapor, momentum and continuity equations were developed for 


the isothermal vapor flow. 





Figure 1. Coordinate system and geometry 
for condensate film analysis. 


dEah 





A. CONDENSATE MOMENTUM EQUATION (X-Direction) 
By assuming that the momentum changes in the condensate 
are small a static force balance may be taken on an infini- 


tesimal liguid element in the x-direction as shown in 


Figure 2. 

0 = st - SB + ppw’ rsin¢ (1) 
where 

tT = shear stress 

p = pressure 

xX = co-ordinate measuring distance along surface 

y = co-ordinate measuring distance normal to surface 

r = radius 

@ = half cone angle 


oa density of liqwid 


W = angular velocity 


B. CONDENSATE MOMENTUM EQUATION (Y¥-Direction) 
Similarly, using Figure 2, a force balance in the 


s 


Y-direction may be taken 
een : 
0 a P-w° x cos-¢ (2) 


Making use of Figure 2, the radius at any point in the fluid 


condensate may be expressed as: 
es) = hae ex San d -y cos $ (3) 


where 


Ry = minimum wall radius in the condenser section 


12 








Figure 2. Cross-section of infinitesimal liquid 
element. 
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SW@estituting equation (3) for r(x,y) into (2), and inte- 


Seating (2) from p to P,, and y to 6 yields: 
p(x,y) = py(x) + p,w°cos ${(R,+ x sind) (8-y) 


= $ Be Gi See (4) 


Differentiating equation (4) with respect to x gives 


dp 
: ae as 
ae a a + p,w'cos o{(R,+ x sin $-écos $) F> 


+ (6-y) sin 6} (5) 


Substituting equation (5) into equation (1) and integrating 


from y to 6 and T to Le gives 


dn aie 
T=- = (6-y) - T- P,w*cos ${(R +x sind-dcos $) S"(6-y)} 
+ pwsing{R, +xsin$-6cos ¢}6-y) 
Re du (x,y) 
where r 


Us Viscosity of the liquid 


u(x,y) = velocity of the liquid 


Integrating equation (6) from 0 to y and from 0 to u gives 


an expression for the velocity of the fluid: 





dp Z Tey, 
: 1 V y V 
TOR A) SS eee [i === 
p gut asl, v7 
+ [ + x sin¢-Scosd|[sing-cosoSs fy ~ S| 
Dele L 4 
(a) 
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This expression for u(x,y) agrees with Ballback's except 


for the shear stress term. 


Ge CONTINUITY EQUATION 
It is apparent from Figure 3 that due to the conserva- 


eon Of mass, with condensation 


dw. dw. 

a ae (10) 
or 

We = Nes . (10a) 
where 

We = Mass flow rate of fluid 

w., = mass flow rate of vapor 
The mass flow rate of the liquid may be written as 

6 
We = : Pu 2mr dy (11) 


The integrating of equation (11) using equations (3) and 


(7) yields 


3 


We = 2719, t- We ax (Rt x sin 9) fon COs > aq 6 , 


le 52 
+ 210 ¢ {- et (Ry Sle ee COo oe 84) 


210 --w~ 
(Ry +xsino- Scos¢(sing - cos ¢ $ oS) 
if 
(R Hxsin OS: = cee ee (12) 
[ a x<Sin 3 Cc 74 


The mass rate of flow of the vapor may also be written 


LS 








av OM yy ay 
ae = ax) w+ aa ax) 
+ aaa ae 
Py ax V 
AC 
V 
ene 

WwW aw 

i a 
ie aS ax 


s section of infinitesimal 


Figure 3. Cros 
nad vapor elements. 


ji quire a 
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i = p.,7R V (3) 


where 


R= Ry +xsin 6 - S6cos ¢ 


D. ENERGY EQUATION IN CONDENSATE 
From Figure 2, using an energy balance, the rate of 


condensation may be written 


dQ _ 2 dw. - 2n(R\+ x sin ¢) (T_-T,) ae 
dx fg dx Oe caer = 

k k h 

f W 
where 

Te = saturation temperature of vapor 
T., = inlet coolant temperature 
t = thickness of condenser wall 
6 = film thickness of condensate 


ke = thermal conductivity of condensate 


Le = thermal conductivity of condenser wall 
h = outside condenser wall heat transfer 
Ceciricient |. 


E. VAPOR MOMENTUM EQUATION 
Using Figure 3, and neglecting the momentum of the 
condensing vapor, a force balance may be taken in the 
x-direction. 
dp 
A oie) (Aco =~) dA - TGA, 


Py 


= (w t+dw_) (vtdv) SEN (15) 


where 


ey 





P.. = pressure of vapor 


Vv 

v = velocity of vapor 

i = Cross sectional area of vapor element 
AY = surface area over which the shear 


stress acts 
Neglecting higher order terms, and using the fact that 
dA = 2TRdxX equation (15) reduces to 
Vv o dv Vv d Ty 


A: ee rae aaa ee “R | (16) 


Using the expression TS : pees where f£f = friction factor, 


and substituting into equation (16) yields 


Z 








dp Oe £ 
Vo gy RL va 2 v 
dx PY ax Bee Gs (ey ee oe) 
Expanding the derivative in equation (17) yields 
2 2 
on a OO ai a i (17a) 
ax Vv -dx R dx R 


The first two terms in equation (17) represent the momentum 
effects of the vapor while’ the last term represents the 
usual Preeti on effect. For the low condensate velocities 
which were expected during the experimental runs, according 
to Bergelin et al. [2] the friction factor is approximately 
equal to the friction factor used in smooth tubes (that is, 
the effect of the liquid film on the condenser surface is 
negligible). According to Rohsenow and Choi [3], the 
Beretiion factor for turbulent flow in a smooth tube for vapor 


Reynolds numbers between 2000 and 50,000 is expressed as: 


18 





20791 AIO I/SEE 


= (Re) CES) 
Ey 


F. DETERMINATION OF FILM THICKNESS AND HEAT TRANSFER RATE 


/ (Re) * 





nauatrons (10), (12), (13), (14), .and (17a) canbe 
further reduced to two simultaneous ordinary first order 
differential equations. Using equation (14) with the use 


dw 
Oimequations (10) and (13) for — yields: 


Aa 2(R, +x sin $) (T, =) 


QV aé 
=a ai aa [sine cos > =| ee 
aie: MR h fo 

W 


At the same time, using equation (12) in conjunction with 


(10a) and (13) we obtain 


a 3 
ae ( 7 Py, . 0 eee ae) ra 
PTR V= TES Te ae +x sing 1 cOs | 74 u ‘| 
1 32 53 
+ 270, fio ee ot cos $ >] 
22 
270 —¢ W , ; ds, 
+ ee (Ro+ x sin$- cos@¢) (sin¢- cos a 
Rapeecis a aoa hy = 8) (20) 
EX a x sind 3 a4 


If equation (17a) is now substituted into equation (20), an 
equation for ce is obtained: 


ax 2) 
Py Po i. dv 
2 a 
pey. bara ap 8) tear = 2p - 
as 2 Vv 20¢ Tax 





dé _ 4 Pee eee ee eS 
P COs Pw 2 
where 
3 
: 6 5 4 
Py = es x sin?) z= = cos $ a7 4 
62 63 
ays (Ro+ x sin 4) = 7 cos $ = 


ie 





Equations (19) and (21) are two simultaneous first 
order ordinary differential equations with unknowns 6 (x) 
and v(x). A Runge-Kutta numerical integration was used to 
solve them with the IBM digital computer. It was hoped 
that by using first order differential equations rather 
than a second order, as Daley did, that the sensitivity to 
initial 6 would be reduced. Figure 4 is a flow chart of 
the program used to solve the differential equations. 


The integration was started by assuming initial con- 


ditions, b=65 and v=0 at x=0. With these initial conditions 
ae = tangd at x=0. Figure 5 shows a schematic representing 


the calculated 6(x) for different initial values of 6. 

Note that if 65 is too small the calculated 6 becomes 
negative before the condenser overfall (physically impossible 
Sorution). If 65 is too large the calculated § grows at a 
constant slope equal to tan$?(physically unrealistic). With 
the proper choice of és the calculated §6 reaches a maximum 
along the length of the condenser and then tapers off to a 
value close to zero at the overfall. The exact value of 6 
at the overfall is not known. It was observed however that 
the calculated value of 6 at the overfall was approximately 
a tena Oo, and the heat transfer rate was insensitive to an 
exact value at the overfall. This result agrees with Nimmo 
and Leppert [8] who experimentally found for laminar film 
condensation on a finite’ horizontal surface that the heat 


transfer rates were practically independent of the end 


@enditions. 
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INPUT DATA 


Tor Por Per Up rl Ke 


w! h, ae Sy5 


t, RPM 





CALCULATION OF Q SUBROUTINE FRI 
Using the value of 6 found in Numerically integrates 
Subroutine RKl equation (14) is equations (19) and (21) 
used to find the heat transfer using 1000 intervals 
rate at each interval of x. Q Cf. rom x=0 ro x=L 





is obtained by using the trap- 


ezgoidal rule. oe 


SUBROUTINE DFERS 


Contains explicit. «une 
tions of derivatives, 


PRINT 6, v, O dv ais 
cx ax” 
Figure 4. Flow chart Of computer program. 
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2 


2 S25 a a 
Vee ee eee oo 
Py E 2 = PYV 20. 2Py | 





tan ¢ = —————______ 
¢ 5 (22) 
5 2p.Vv 
Py COS $)/P -w Re 
2(R_+ - 
dv _ ( o" x sind) (T. ie) _ 2vsind oo 
ax ae See! Z R 
Pale | Se, 
W ; 


Equation (22) may be further simplified by an order of mag- 
nitude analysis. Neglecting the second terms of Py and P5 


and using the expression, R = Ry tex Sind With) tite 

UeRZ Pp p 
dv £ 1 2 
So << ae oS 
ax 20 - ge! poe 


equation (22) may be rearranged in terms of a cubic equation 


previously observed result that Py 


tore O (XxX) : 





oe 
5 VATOy a NA Z 
oe 2 V - 2 6 
3 R sind Pew R + = pV a ro 
Z 
UeR 





= 0 (24) 
18 


Equations (23) and (24) were solved simultaneously for 
v(x) and 6(x). Equation (23) was numerically integrated 
using 1000 steps of integration over the cooling length of 
the condenser. The integration was started with initial 
conditions és = 0 and v=0 at x=0. Having found v(x) at a 
particular step, equation (24) was solved by a subroutine. 
This value of 6 was used to find v(x) for the next step of 
integration. The results of approximation 1 are recorded 
tmerapile £ for half cone angles Of 131° ;,.2.2°, ends.3. 

Using approximation 2 on equation (21), equation (24) 
was further simplified and rearranged to yield the. followi'g 


expression for 6(x): 





TABLE I 


DHEOHEAT FLOW, QO (Btu7hr), FOR 
DIFFERENT HALF CONE ANGLES 


NUMERICAL APPROXIMATE APPROXIMATE 
> SOLUTION SOLUTION 1 SOLUTION 2 
Eqs. (19) & (21) Eqs. (23) & (24) Eqs. (23) & (25) 





0° 8,505 = “ 
ae 10,082 9,836 10,290 
2° LE aie 11,033 ak 292 
° i 910.4 11,769 11,998 
3P UV WES 
(68) oe we (25) 


20. Hsin o 


Using equations (23) and (25) a numerical solution was found 
for v(x) and 6(x) in the same manner as the solution for 
approximation 1. The results are tabulated in Table Teen 
feweeeconce angles of .1°, .2°, and .3°: 

Table I shows that approximation 1 is lower than 2. 
This is to be expected since approximation l includes both 
interfacial shear between the’ vapor and condensate and the 
pressure drop in the vapor. At the same time approximation 
1 is less than the more accurate solution of equations (19) 


enamel), thus giving a conservative result with a deviation 





*tn finding the solution for .3°, 10,000 steps were 
used in the integration. This was necessary because as the 
half cone angle became larger than ,2°, the numerical ime. 
gration became very sensitive to the step size. 
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of only 2.5%. Conseguently, approximation 1 was used to 
arrive at a theoretical solution for 6(x) and Q for the 
one degree half cone angle condenser used in the experi- 
mental program. The results are plotted in Figures 9 


ana 10. 
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II. EXPERIMENTAL PROGRAM 


A. DESCRIPTION OF EQUIPMENT 

In order to facilitate the description of the entire 
experimental system, the system is sub-grouped into four 
parts: (1) evaporator, (2) condenser, (3) instrumentation, 
and (4) auxiliary equipment. Included in the description 
are drawings which are intended to aid the reader in 
visualizing the different components of the system. 

Figure 6 is a schematic representation of the entire system 
including the auxiliary equipment. Figure 7 is a cross- 
sectional drawing of the assembled heat pipe. Figure 8 
is a photograph of the heat pipe and associated equipment. 
lL. Evaporator 

The evaporator, modeled to approximate the rotating 
boiler tested at Lewis Research Center [17], is made of 304 
stainless steel. It has an inside diameter of 3.125 inches, 
an outside diameter of 3.635 inches, an exit diameter of 
Peo25 inches, and an overall length of 5.905 inches: 

The actual heated length is 3-7/8 inches around which 
an ll-gauge Chromel-A heater wire in a 3/16 inch inconel 
sheath is helically wound. The procedure of bonding the 
heater wire to the evaporator is an exacting process in 
order to ensure a uniform heat flux axially along the heater 
surface and to guard against electrical shorting between 
the heater wire and other components. The procedure was as 


EoUMOWwS : 
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a. A sheet of gold-nickel foil, .003 inch thick, 
was tack welded to the outside evaporator surface. 

b. The heater wire was helically wound over the 
gold-nickel sheet with two wraps of .020 inch diameter Golde 
nickel wire wound intermittently between the heater coils. 

c. The assembly was then brazed in a hydrogen 
eaumospherically-controlled oven. 

d. Sauereisen cement was baked around the outside 
of the heater wire using the heater wire as the heat source. 

e. Asbestos insulation was wrapped around the 
heater wire and tied off with glass tape. 

f. Fish paper was placed into the heat sink on the 
condenser end of the evaporator to insulate the heater 


coils from the evaporator. 


Power was furnished to the heater coils by a DC motor- 
generator capable of delivering 150 amperes at 220 volts. 
The electrical power was passed hy a brush assembly con- 
Sisting of eight spring-loaded graphite brushes through two 
beonze collector rings to She helically wound heater wire. 
The rings were supported by twelve 1/4 inch stud holts. 
Phenolic spacers were used to align the collector rings and 
to insulate the bolts from the rings. 

In order to reduce the axial heat losses, two radial 
1/16 inch wide grooves were machined in the outside surface 
to within 1/16 inch of the inside boiler wall at either end 


of the heating element. 


30 





The outside face of the evaporator is closed by the 
insertion of two pyrex glass view Ports, 3=l72™inenes in 
diameter. The inner view port rests on a neoprene Shovel yeXs; 
seal and is separated from the outer view port by a 
3/32 inch compressed fiber gasket. The view ports are held 
in place by a stainless steel end flange secured in place 
by twelve 1/4 inch stud bolts. The original design of two 
solid glass view ports was changed hy fusing the inner view 
port to a glass tube which passes through the outer view 
port. This tube allows a vacuum to he drawn from the inside 
of the heat pipe during the filling procedure, and is easily 
sealed. 

The inside diameter of the evaporator is stepped 1/2 
inch from the exit diameter to allow for the formation of a 
liquid annulus in the evaporator. 

Ze Condenser 

The condenser section is machined of stainless steel 
in the shape of a truncated cone. The half cone angle is 
1°. The base of the condenser mates to the evaporator exit 
and has an inside diameter of 1.809 inches and tapers over 
a 10 inch length to a diameter of 1.460 inches at the 
truncated end. At this point, the condenser extends an 
additional 1-1/2 inches to house a cylindrical end plug. 

A step is machined at the evaporator-condenser inter- 
section where the condensate flows from the condenser to the 
recessed inside diameter of the evaporator. The first inch 


of the condenser surface near the evaporator is assumed to 
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be adiabatic and 9 inches of the condenser's overall Fengith 
is considered as the cooling surface. The wall thickness 

of the condenser is 1/16 inch in order to reduce its thermal 
resistance. 

The condenser section has two flanges. The large flange 
mates to the bearing seat flange on the evaporator and is 
held in place by the 1/4 inch stud bolts. The other flange 
Mates to a flange on the end plug held together by eight 
1/8 inch stud holts. The end plug has an extension on its 
open end on which are mounted four Phenolic juNeELOnN boards. 

3. Instrumentation 

A schematic drawing of the instrumentation set-up 
Re shown in Figure 6. 
a. Temperature Measurements. 

Temperature measurements are obtained from nine 
grounded copper-constantan thermocouples and two quartz 
thermometers. The thermometers are used to Measure a 
temperature profile in the evaporator wall, the liquid 
annulus temperature in the evaporator, the vapor saturation 
temperature in the evaporator, the outside wall temperature 
of the condenser, and the temperature at the location of the 
pressure transducer. 

In order to measure the temperature gradient in the 
evaporator wall, four thermocouples with a 1/16 inch dia- 
meter stainless steel sheath are utilized. Fach thermo- 
€oupte is placed within a 1/16 inch well drifted 3-872 


inches deep into the evaporator wall. The thermocouple 
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junction is coated with epoxy to give it enough strength so 
that it will not bend when placed against the end of the 
well. The epoxy is sanded such that the junction makes 
contact with the wall. 

One thermocouple is located 1/16 inch radially from the 
inside of the evaporator, another is located 1/16 inch 
radially in from the outside evaporator surface, and the 
remaining two are located opposite one another in the middle 
of the 1/4 inch wall at the same radial position to check 
on the uniformity of the temperature field. 

The leads from all of the thermocouples in the evapor- 
ator wall are brought through slots in the evaporator and 
forward condenser flanges, down the side of the condenser, 
and through the rear condenser flange to phenolic junction 
boards. From the boards, the leads go inside the shaft and 
down to the slip ring unit which is attached to the main 
assembly by a coupling mounted at the end of the shaft. 

A thermocouple is suspended in the liquid annulus to 
measure the liquid bath temperature. The remaining thermo- 
couple in the evaporator is suspended in the vapor space to 
measure the saturation temperature. 

Since these two thermocouples protrude into the inside 
of the heat pipe where a vacuum is held, precautions were 
taken to ensure that a seal was maintained. The ends of 
the two thermocouples were enclosed in stainless steel 
flanges and the leads were soldered at the slots on the 
condenser forward flange to keep air from passing through 


the slots into the evaporator. 
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In order to measure the outside condenser wall temper- 
ature, two thermocouples were soldered to grooves on the 
outside of the condenser at locations of 1/3 and 2/3 of its 
length respectively. 

The remaining thermocouple was used to measure the 
temperature at the pressure transducer. The junction is 
soldered to the inside of the pressure transducer arm 1/32 
inch from the surface. 

The temperature of the coolant is measured before and 
after spraving by quartz thermometers mounted in mixing 
tubes located above and below the condenser spray: cooling 
cylinder. 

b. Pressure Measurements 
The saturation pressure of the vapor at the 
small end of the condenser is measured by a semi-conductor 
pressure transducer mounted in the vapor space. The trans- 
ducer is temperature compensated and is threaded on a 1/4 
inch diameter arm that extends from the shaft through the 
end plug and about three inches into the vapor space. The 
leads of the probe pass through the transducer arm up 
through the shaft to the junction boards. From there, the 
leads pass back through the shaft to the slip ring unit. 
c. Coolant Flow Rate Measurements 
The coolant flow rate is measured hy passing 
the coolant through a flow rotameter before the coolant 


enters the spray nozzles. 
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d. Rotational Speed Measurements 
The rotational speed of the heat pipe is 
measured by two meters. A Beckman Eput Meter measures the 
RPM using a magnetic pickup, and an electronic strobe scope 
is used as a check. 
e. Electrical Power Measurements 
Power to the heater coils is measured by an 
ammeter and voltmeter connected in series and parallel 
respectively with the power lines returning from the heater. 
The power level is controlled by a field-rheostat placed in 
parallel with the DC motor-generator. A 150 ampere circuit 
breaker is included in the circuit. 
foo Ring Unit 
The slip ring unit allows the rotating thermo- 
couple leads and pressure transducer leads to pass electric 
current to stationary leads. The slip ring unit used in 
the experiment utilizes liquid mercury as a viscous contact 
medium between the rotating and stationary leads. The 
frictional heating noise is less than lO microvolts at 
speeds up to 4000 RPM. The slip ring unit has 22 terminals 
to accommodate the nine thermocouples and pressure trans- 
ducer. Sthe output Exveom the slip rings 1s fed to a Hewlett- 
Packard 2010C Data Acquisition System containing an inte- 
grating digital voltmeter, guarded data amplifier, and 
guarded cross-bar scanner. This system provides printed or 


Vike moutputL with anvwaccuracy of + .5 microvolts. 
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4. Auxiliary Fquipment 


The auxiliary equipment includes: (a) drive 
assembly, (b) test stand, (c) safety shields, and (da) spray 
assembly. 

a. Drive Assembly 

The drive assembly consists of a drive motor, 
end plug, shaft, pulley and bearing. The drive motor drives 
the shaft by means of a 2.65 inch diameter pulley and 
V-belt. The driving motor is a 2 horsepower, 3-phase 
variable speed motor capable of a maximum speed of 4000 RPM, 
The speed is controlled by a pushhutton electric remote 
control unit and the motor can be shut down instantaneously 
by a magnetic disc brake. 

The end plug has a 3/4 inch threaded shaftway machined 
into it. Between the end plug and the shaftway is a .27 
inch hole to allow for the insertion of the pressure 
transducer. The stainless steel shaft is hollow to allow 
for the passage of the thermocouple and, the transducer leads 
tortne slip cing unit. The shaft has 3/4™imeh Of “rts 
outside diameter threaded for screwing into the shaftway, 
and is threaded for 3/8 inch along its inside diameter on 
both ends. The condenser end mates with the transducer arm 
described earlier and the opposite end is joined to a 
coupling connected to the slip ring unit. 

There are two main hearings supporting the heat pipe. 
Mie large, 5.96 inch diameter, floating, single bearing 


seats on the evaporator flange which is supported by an 
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aluminum pillow block. Copper cooling lines through which 
tap water flows have been placed in the pillow block to 
help cool the bearing. There is an oil input line at the 
top of the pillow block which allows for lubrication of the 
bearing. The small double row, angular contact hearing on 
the condenser end is mounted on a seat machined on the shaft 
and is self-lubricated. 
b. Test Stand 

The test stand was designed and built so that 
the heat pipe could be rigidly supported with its longi- 
tudinal axis ranging in position from 0 to 90 degrees from 
the horizontal. Both the heat pipe and variable drive 
motor are mounted on steel hbed-plates connected by a 4 inch 
iron pipe. The iron pipe is held in place by three 2-inch 
thick clamps which are supported two feet off the ground 
by a 1/4 inch steel plate. As the heat pipe is rotated in 
any position, the drive motor rotates with it and the entire 
assembly is then securely held in place. The motor has heen 
vibration mounted in order to keep the vibration to a 
minimum. 

c. Spray Cooling Assembly 

Cooling of the condenser section is done hy 
spraying tap water on the outside of the condenser during 
rotation. The condenser section is enclosed by a 9-inch 
long 13-inch diameter stainless steel cylinder, the con- 
denser box. The ends of the condenser hox have Garlock 


rubber seals separating them from the condenser section. 
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The condenser box is slid into place over the condenser 
section through its open ends and is then bolted in place 
to the support plate. An 8-inch x 6-inch window has been 
cut in the condenser box in order to be able to make 
adjustments or repairs without removing the condenser box. 
During operation, a steel plate is holted over the window. 
The condenser box has eight spray nozzles mounted at the 
same axial position. The nozzles may he varied to give a 
range of droplet sizes from 300 to 600 microns. 

The coolant flows from copper feed lines to a mixing 
tube located at.the top of the concenser box. A quartz 
thermometer records the temperature of the coolant in the 
cane, From the tube, the coolant flows through plastic 
tubing to each of the spray nozzles. The coolant drains 
from the bottom of the condenser box to another mixing tube 
where the temperature is recorded hy the other quartz 
thermometer. From there the coolant flows to an external 
relaarc aly alge 

ad.) Ssatety Shields 
During operation, the entire heat pipe assembly 
is surrounded by 1/8 inch thick stainless steel shielding 
bolted to the support plate to ensure safe operation. 

The evaporator shield, shaped like an inverted U, is 
bolted to a 1.0 inch thick plexi-glass window on its open 
end. This window serves as a safety shield and also allows 
the internal pipe operation to be observed. The other end 


Mates with the condenser box. 
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B. EXPERIMENTAL PROCEDURES 

The experimental procedures consist of: Chee ling 
Procedure and (2) Operational Procedure. 

i eeP liane eproceeure 

In order to obtain film condensation in the conden- 

ser, the inside surface of the condenser must be prepared 
so that proper "wetting" will take place during operation. 
The following procedure is used: 


a. Scrub the condenser surface with’ 190°P Albeconox 
detergent solution. 


b. Flush with distilled water. 

c. Spray with alcohol. 

a. Rinse with distilled water. 

e. Spray with acetone. 

£. Rinse with distilled water. 

Having completed the above procedure, the heat pipe is 
ready to be filled. In this study, distilled water was used 
ao othe working fluid. 

In order to get a sufficient range of Saturation tem- 
peratures and not exceed the 35 psia limit on the heat Pl per 
distilled water ice was placed in the evaporator and a low 
vacuum was drawn. The filling procedure was as follows: 

a. System was cooled by circulating cooling water 
through the condenser nozzles and by inserting distilled 
Water ce in the evaporator. 

b. Six distilled water ice cubes (approximately 
150 ml) were subcooled hy flashing them With LouLd Ny and 


placed into the evaporator With tongs. 
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c. Glass windows were put into place and the flange 
was tightened with a 30 in-lb. torque. 

d. The vacuum was drawn through the glass tube 
connected to a mechanical vacuum pump. 

e. While the vacuum was being pumped, the glass 
tube was sealed using a gas flame. 

2. Operational Procedure 

Having completed the filling procedure, the heat 
pipe is ready to run at a pre-determined RPM. About 10 
a@zops of lubricating oil are added to the oil fill tube 
at the top of the large pillow block. The cooling water 
lines are turned on so that the rotameter is at approxi- 
mately 20% flow. 

The drive motor is turned on and the heat pipe is 
slowly brought up to the desired RPM. Once a liquid 
annulus has formed in the evaporator and the desired RPM 
has been reached, the DC motor generator is turned on and 
power is added to the evaporator. Having reached thermal 
equilibrium for a particular power setting, data was taken 
by the instrumentation described earlier. 

Shut-down is accomplished by first shutting off the 
current to the heater, and then decreasing the RPM by use 


Olmunen remote Control Unit. 


Coe EXPERIMENTAL RESULTS 
The heat pipe was run at 702 RPM and 1404 RPM. The run 


of 702 RPM was chosen because it was the lowest RPM at which 
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centrifugal forces dominated the gravitational effects. 
Power inputs to the heat pipe ranged from 1 kW to 9 kW. 
1. Preliminary Results 
Using distilled water as the working fluid, the heat 
pipe was rotated at 702 RPM to obtain the heat flow through 
the evaporator wall for a particuler power setting. The 


method used to compute Q 


ae the heat flow through the 


evaporator wall, and the power supplied to the heater coils 
is shown in Appendix B. Table II shows a comparison between 
the electrical power into the heater coils and the heat flow 
through the evaporator wall for that particular power 
setting. 

The data taken from the four thermocouples in the 
evaporator wall is obviously incorrect since part of the 
data shows that the heat flow through the evaporator wall 
is higher than the electrical heat in for that particular 
power setting. In other cases it appears to he low. 

There are two explanations for the incorrect measure- 
ment of the heat flow through thesevaporator wall. The 
first is the fact that the reading of thermocouple #2 
appeared sei When the heat pipe is NOEmroOtLatand , 
thermocouple #2 appears to be reading correctly, but when 
the heat pipe rotates it begins to read inadequately. It 
was concluded that the problem existed in the slip ring and 
may be due to insufficient mercury in the slip ring. The 
second explanation pertains to the reason why the heat flow 


through the evaporator wall is higher than the heat into 
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TABLE II 


POWER IN AND HEAT FLOW THROUGH FVAPORATOR WALL 


HEAT FLOW IN 


eee, EVAPORATOR WALL 

(BTU/HP) 

3,345 4,059 
6,712 10,315 
9,963 12,848 
13,325 7,610 
16,988 9,126 
20,751 24,041 
24,031 13,437 
D7 278 31,986 


errr erm mm mn 


the heater coils for a particular power setting. It was 
concluded that the outside edge coils on each side of the 
heating surface were conducting axially. As a consequence, 
there wasn't a uniform temperature distrihution axially 
along the heating surface rae assumed. Therefore, the actual 
heating surface length was less than the assumed length in 
the calculation of the heat flow and would account for the 
high values of the heat flow obtained in Table II. 

At this point, it was decided to obtain the heat flow 
out of the heat pipe by performing an energy balance on the 
condenser cooling section. The coolant flow rate was 


determined by a flow rotameter, and its mean temperature was 


Meas@Ecd before and after spraying, by quartz thermometers 
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mounted in mixing tubes located above and helow the conden- 
ser box. Having this information, along with the coolant's 
specific heat, the heat flow was determined hy a simple 
energy balance. In order to make this method more accurate, 
the inside of the condenser box was insulated with silicone 
foam rubber and the plastic coolant lines were insulated 
with fiberglass. Appendix B contains a table showing the 
heat flow out of the condenser at different values of the 
heat pipe saturation temperature. 

2. Visual Observations 

Using a strobe light visual observations were made 

through the view ports. At 702 RPM, the light provided by 
the strobe light was insufficient. To increase the light, 
the strobe light's frequency was increased to three times 
the RPM of the heat pipe. As a consequence, the ohservations 
made at 702 RPM were not clear. The obServation of the 
liquid vapor interface in the evaporator at low saturation 
pressures indicated violent boiling with large vapor 
bubbles breaking through the Pequidyvaporsinterface at 
702 RPM. At higher saturation pressures, the boiling was 
more controlled with small bubbles constrained to the 
surface of the liquid. At the same saturation pressure, the 
boiling was less violent at 1404 RPM. 

The visual observation in the condenser indicated that 
there were ripples and random hreak-up of the condensate 
film, exposing parts of the condenser surface to the vapor. 


Furthermore, the condensate swirled and washed the condenser 
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wall at a higher RPM than the wall itself. It is doubtful 
that the above observed effects were due to condensate 
turbulence. According to Rohsenow and Choi [3], in the 
presence of interfacial shear, turbulent condensate flow 
occurs for a Reynolds number between 100 and 1800. For the 
low vapor velocities which occurred during the observed 
experimental runs, the calculated transition Reynolds number 
in the condensate was approximately 1000. However, the 
highest Re calculated for the condensate was only 291. It 
was concluded, therefore, that the swirling and random film 
break-up were not caused by condensate turbulence hut were 
probably caused hy vapor interaction with the condensate. 

| Except for the 1 kW case, the Reynolds number calculated 
for the vapor flow for all of the power inputs was hetween 
2000 and 9000 (for 1 kW, Re = 972). Since flow in a smooth 
tube is turbulent for Reynolds numbers between 2000 and 
50,000 [3], it was concluded that the vapor flow was tur- 
bulent except for the 1 kW case. The swirling of the 
condensate: was probably caused by the higher velocity vapor 


striking the condensate. 
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Pil. SDISCUSSIONSOFR RFSULEIS 


The graphs in Figures 9 and 10 show the experimental 
and theoretical heat transfer rates as a function of the 
saturation temperatures (corresponding to an internal vapor 
pressure). 

In order to find the theoretical values of the heat 
transfer rate, the value h, the heat transfer coefficient 
was needed. The method used in determining h is shown in 
Appendix B. Because of experimental uncertainties, the 
values of ayy = 1000 and eae = 1800 were used since these 
values were very close to the average of the minimum and 
maximum experimental values of h determined for the dif- 
ferent operating conditions. 

As observed in Figure 9, the theoretical values of On 
and Oe. are less than the experimental value of Q at 
702 RPM. In the development of the analytical solution it 
was assumed that smooth film condensation occurred in the 
condenser. The visual observations described in the experi-~ 
mental results indicated that a smooth condensate fiim did 
not exist for the operating conditions used in these 
experiments but rather that the condensate film swirled and 
washed the condenser wall and appeared to hreak up in a 
random manner expoSing part of the condenser surface to the 


vapor. The swirling effect and the resulting increased 


Gomvecerton in the condensate Increased the heat transfer 
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Figure 9. Heat transfer rate vs saturation temperature 
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rate. It was also possible that the ripples and film 
break-up decreased the overall film thickness, reducing the 
thermal resistance, and thereby increasing the heat transfer 
Eace. . 

Figure 10 shows that the analytical and experimental 
results compare more favorably for the 1404 RPM run. There 
are two probable explanations. First of all, the higher 
RPM could have flattened the film causing it to approach a 
smoother film than in the 702 RPM run. If this were the 
case, the smooth film condensation assumptions would be 
better, and the.analytical solution would he closer to the 
experimental results. Although the visual observations 
‘seemed to indicate that there were no non-condensibles 
present, the presence of non-condensibles in the condenser 
could be the second explanation. If non-condensibles were 
present, the experimental heat transfer curve would he 
shifted to the right since the heat transfer rate would 
require higher vapor pressures (and, therefore, higher 
saturation temperatures). — | 

Lf the inlet coolant temperature were held constant, 
the heat Gee ee rate should increase with increased RPM 
for the same saturation temperature. For example, the 
analytical solution for oon at 1404 RPM at approximately 
ete Ssameet. used at /027ReM resulted gaa higher heat transfer 
rates as shown by the dashed plot in Figure 10. Note that 
the experimental data at 1404 RPM does not appear to show 


an ine@rease in heat transfer rates from the 702 RPM run. 
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This is due to the fact that the 1404 RPM run was conducted 
at a hagher T. (7° difference) than the 702 RPM run. This 
increase in T, increased the saturation temperature by 
approximately the same amount such that the experimental 
heat transfer rate data points were shifted to the siighis. 
At the same saturation temperature, therefore, the heat 
transfer rate was lower for the higher T, even though nate 


was at higher RPM. 
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Iv. CONCLUSIONS AND RFCOMMFNDATIONS 


A. CONCLUSIONS 

1. The rotating non-capillary heat pipe is an effective 
heat transfer device. 

2. “Smooth film condensation was not observed during the 
experimental runs. 

3. The approximate numerical solution conservatively 
predicts the heat transfer rate with a deviation of 18% at 
702 RPM. ‘The deviation of the numerical SOlUtION decreases 
to 5% for 1404 RPM. 
| 4. The performance of the heat pipe is strongly 


dependent on the method used to cool the condenser. 


B. RECOMMENDATIONS 

1. Run the heat pipe at higher RPM to determine 
whether increased RPM provides better agreement between the 
analytical and experimental results. 

2. Increase the instrumentation of the condenser to 
give more accurate data for the experimental heat transfer 
rate out of the condenser. 

3. Revise the analytical solution to include the ripples 
and film break-up in the pene enee 

4. Obtain experimental results for the effect of non- 
Pesndeneibles in the heat pipe. 

5. Use different methods of cooling the condenser wall 


by using different spray nozzles. 





6. Use a condenser section with a half-cone angle 
greater than 1° to obtain experimental data for the effect 
Seelarger cone angles on the heat transfer rate. 

7. Use working fluids with boiling points lower than 
water in order to find optimum heat transfer rates at lower 


saturation temperatures. 





APPENDIX A 
CALIBRATION 


1. Calibration of Pressure Transducer. The pressure trans- 
ducer consists of a temperature-compensated, senticonductor 
strain gage Wheatstone bridge mounted on a 1/4-inch diameter 
diaphragm. It was calibrated against a Wallace-Tiernan 
standard pressure gauge (with a known accuracy of 0.03 + 
0.01 psia) for pressures greater than atmospheric, and 
against a mercury manometer (with eimmeectainty Of et .02 
psia) for pressures less than atmospheric. The transducer 
was calibrated at different operating temperatures by 
placing it within a controlled temperature oven with a 
quartz crystal thermometer (with an uncertainty better than 
20 Lo E). Reproducible data was obtained from 0-30 psia 

and from 65-222°F, and the transducer sensitivity was 
measunwed to be 3.54 + 0.06 mv/psia. Although the transducer 
Wheatstone bridge was temperature compensated, a non-linear 
zero shift with temperature was measured and incorporated 


eo a Sserres Of calibration curves. 


2. Calibration of Thermocouples. Nine 30 gauge, Kapton 
insulated, copper-constantan thermocouples were carefully 
manufactured at the Naval Postgraduate School. They were 
calibrated up to 384°F using the same Penper Seuie controlled 


Oven and quartz crystal thermometer mentioned ahove. In 
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addition, the following fixed point calibrations were made: 

a.) Gripile point of water, (32-0 F) 

b. Boiling point of water (212.0°F at 14.696 psia) 

c. Melting point of tin (449.4°F) 

d. Melting point of lead (621.3°F) 

There were a total of nine calibration data points 
taken over the range of 32-621.3°F. Using the NBS table 
of calibration data for copper-constantan thermocouples as 
a zero reference, a deviation curve of the error (mv) vs. 
the thermocouple reading (mv) was plotted. This curve (with 
an uncertainty of 0.1%) was then used as a correction 
factor for each data point taken during the experiment 
‘before entering the NBS table for the temperature. Recali- 
bration data was taken at a later time and the repeatability 
of all nine thermocouples was excellent. Using the boiling 
point of water and the melting point of tin, the deviation 


was within 0.005 mv of the previous data points. 


3. Calibration of Rotameter. The method used in the cali- 
bration Be ne flow rotameter consisted of using a Toledo 
scale to measure the mass of the coolant flowing into a 
bucket placed on the scale. As coolant flowed into the 
bucket, an arbitrary starting point on the scale was taken 
and a stop watch was used to Wane the time it took for 
the weight to reach 20 pounds. This procedure was repeated 
several times for each per cent flow rate used on the 
rotameter. Data points were obtained for 14, 16, 18 and 
19.7 per cent flow rates. From this data a llanear calibra 


tion curve was obtained. 
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APPENDIX 5S 


DATA RFDUCTION 


ie = Calculation, of Qevap 


The heat flow through the evaporator wall may be reduced 
to a one dimensional steady state problem if it is assumed 
that there is no axial heat conduction and steady state 
conditions exist. Under these conditions, the well known 
conduction equation in cylindrical coordinates may he 


weitten [1]: 








la ar| _ 
ee ace 
or 
ar .. ve 
k(T) r a= = B = Constant (a) 
where: 
he (Ob) Ky + m(T- T.) eu 


A plot of k vs. T was made based on data taken from the 
Metals Handbook [5] for the thermal conductivity at different 
temperatures. From this plot the slope, m, and k. were 
determined. 

SUDSEI LU me equation (b) into equation (a) and inte- 


grating yields: 


Eek) = E + m 
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Using a sub-routine for the method of least squares, 
the IBM 360 computer was used to find the slope, B, of 


feee) vs. Inr. From Fourier’s law: 


ae at 
EVAP = L an ke 1G =: | 
or 
Tween le © 
where 
L = the effective heating surface length in the 
evaporator 
° a te) 
EXAMPLE: Toam = 153.7°F, 702 RPM 
Spee a ee" © 
ii — . 320m. ) 
B = 4430. (BTU/hr- ft.) 


Snvap = 9,126 (BTU/hr) 


2. Calculation of POWER IN 
le = Voltage x Current 


IN 


3. Calculation of Fxperimental Q 


OUT 
Sour = ™ C,(T) ~T)) 
where 
m = mass flow rate (lbm/hr) 
C5 = specific heat (BTU/lbm °F) 
Ty = inlet coolant temperature (°F) 
T. = exit coolant temperature (F) 


Mur = heat flow (BTU/hr) 


ote 





PeeOMPERs TSR i=2153.722F, 102 RPM 


Cc. = 1.0 

P 

(T, - To) = G04" + 0. 3) 
Qour = 12,784 + 541 


fee Calculation os hn, OUES nee condenser wall heat transfer 
CoOcrricient 


Sour = BA, IT, ~ Te) 


A= Qour/Ag Ty Te) 


where 
ae = cooling surface area (£t7) 
ibe = outside wall temperature of condenser (°F) 
T., = coolant inlet temperature (°F) 
h = heat transfer coefficient (BTU/hr-ft*-°F) 


EXAMPLE: TSAT = 153.72°, 702 RPM 


Que, = 12,784 4 54% 
AS oe ; 

Cc 

(T - T,) = 28.44 + 6.7 
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5. Tables of Reduced Data 






SAT m 
(°F) 1 (LBM/HAR) 
PSel 635223 
mO7.3)| 832223 
22.8) 804220 
mo7.ot 623224 
moaoeo) 797218 
my2.1)| 764219 
188.3] 764419 
mO4.4) S6321L9 
98.41 870223 
mi5.3) 85024 
mal .8| 825220 
m74.3) 790220 
205.6| Fe6t27 
224.3| 1862427 
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